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Abstract. The strong waveguide dispersion in photonic crystal fibres provides 
unique opportunities for nonlinear optics with a zero-dispersion wavelength Ao far 
below the limit of ~ 1.3 fim set by the material dispersion of silica. By tuning 
the air-hole diameter d, the pitch A, and the number of rings of air holes N, 
the strong waveguide dispersion can in principle be used to extend Ao well into 
the visible, albeit to some extend at the cost of multimodc operation. We study 
in detail the interplay of the zero-dispersion wavelength, the cut-off wavelength 
A c , and the leakage loss in the parameter space spanned by d, A, and N. As a 
particular result we identify values of d (~ 500 nm) and A (~ 700 nm) which 
facilitate the shortest possible zero-dispersion wavelength (~ 700 nm) while the 
fibre is still single-mode for longer wavelengths. 
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1. Introduction 

Photonic crystal fibres (PCFs) P3 E] have led to an enormous renewed interest in 
nonlinear fibre optics pHEJElEI- In particular, the strong anomalous group- velocity 
dispersion has facilitated visible super-continuum generation [7j, which apart from 
being a fascinating and rich non-linear phenomenon jS] also has promising applications 
including frequency-comb metrology [HI El an d optical coherence tomography [TT], 

The regular triangular arrangement of sub-micron air-holes running along the 
full length of the fibre ^1 EI], see Fig. ^ is a key concept for the realization 
of strong anomalous chromatic dispersion which arises in a competition between 
material dispersion and wave-guide dispersion originating from the strong transverse 
confinement of light. At the same time the strong transverse confinement of light 
also serves to dramatically decrease the effective area and increase the non-linear 
coefficient ^3] so that very high optical intensities may be achieved relative to the 
input power. 

The zero-dispersion wavelength may be tuned down to the visible |151 by 
carefully increasing the normalized air-hole diameter d/A while at the same time 
decreasing the pitch A. However, in many cases the very short zero-dispersion 
wavelength is achieved at the cost of multi-mode operation. Furthermore, the 
reduction of the pitch may require a large number N of rings of air holes in order 
to circumvent leakage loss. 

In this paper we in detail study the complicated interplay of the zero-dispersion 
wavelength Ao, the cut-off wavelength A c , and the leakage loss in the parameter space 
spanned by d, A, and N . As a particular result we identify values of d and A which 
facilitate the shortest possible zero-dispersion wavelength while the fibre is still single- 
mode for longer wavelengths, i.e. A c < Aq. 




Figure 1. Cross section of a photonic crystal fibre with N = 4 rings of air holes 
in the photonic crystal cladding surrounding the core defect which is formed by 
the omission of a single air hole in the otherwise periodic structure. The photonic 
crystal cladding comprises air holes of diameter d arranged in a triangular array 
of pitch A. 
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2. Competition between material and wave-guide dispersion 



The linear dynamics and response of a waveguide is typically studied within the 



framework of temporal harmonic modes, E(f,t) = E(f±)e %ut 
where E(r) is a solution to the vectorial wave equation 



V x V x EVr) = e(f,oj)^rE(r). 



(1) 



Here, e(r,u>) = n 2 (r, uj) is the spatially dependent dielectric function of the composite 
air-silica dielectric medium, see Fig. ^ and n(f,u>) is the corresponding refractive 
index. 

Solving the wave equation, Eq. lfT)l. provides us with the dispersion relation 
lo(P) which contains all information on the spatial-temporal linear dynamics of wave 
packets and it is furthermore important input for studies of non-linear dynamics. 
However, quite often the dynamics and evolution of pulses are quantified by the derived 
dispersion parameters /3„ = (d n uj/d n f3)~ 1 with (3\ = v^ 1 being recognized as the 
inverse group velocity. The group-velocity dispersion /3 2 is in fiber optics commonly 
quantified by the dispersion parameter D (typically in units of ps/nm/km) which can 
be written in a variety of ways including 



D 



d 2 p 



dX 



-(3 2 



(2) 



dXduj d\ 2-kc 
where A = 2-kc/lo is the free-space wavelength. The dispersion in the group velocity 
makes different components in a pulse propagate at different speeds and thus pulses 
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Figure 2. Left axis: Dispersion parameter versus wavelength for a photonic 
crystal fibre with d/A = 0.7, A = 700 nm, and TV = 10. The solid curve shows the 
total dispersion, Eq. J5J, obtained through a self-consistent numerical solution 
of the wave equation, Eq. 0, while the dashed line shows the approximate 
addition of waveguide and material dispersion, Eq. 0. The dotted curve shows 
the pure waveguide contribution, Eq. JSJ, while the dot-dashed curve shows the 
silica material dispersion, Eq. JHJ. Right axis: Leakage loss versus wavelength for 
the fundamental and the first high-order mode. The rapid increase in the leakage 
loss of the high-order mode marks the cut-off wavelength. 
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Figure 3. Zero-dispersion wavelength versus fibre diameter for a strand of silica 
in air. The dashed line indicates the cut-off wavelength and the blue shading 
indicates the multi-mode phase. 



will compress or broaden in time depending on the sign of D. The zero-dispersion 
wavelength Ao, defined by D(Xq) = 0, is thus of particular relevance to pulse dynamics 
in general and non-linear super-continuum generation in particular. By pumping 
the fiber with high-intensive ultra-short nano and femto-second pulses near Ao the 
pump pulses will loosely speaking maintain their high intensity for a longer time (or 
propagation distance) thus allowing for pronounced non-linear interactions. 

The finite group-velocity dispersion is a consequence of both the dispersive 
properties of the host material itself as well as the strong transverse localization of 
the light caused by the pronounced spatial variations in the dielectric function. Since 
the variation with frequency of the refractive index of silica is modest (at least in the 
transparent part of the spectrum) the dielectric function satisfies 

e(fx,oj) = e (rx,w) + &(rx,w), (3) 
Se(r±,uj) = e(r ± ,u;) - e (r±,ij) < £ (fj_,w) (4) 

where uj is some arbitrary, but fixed frequency where consequently Se = 0. The high 
attention to the telecommunication band has often made A = 1550 nm a typical choice, 
though this is by no means a unique choice. From e one may define a pure waveguide 
contribution D w to the dispersion parameter by the definition 

where (3q{oj) is the solution to the wave equation with e(f±,uj) — eo{r±,ijj), i.e. the 
frequency dependence of the dielectric function is ignored. Furthermore, this has the 
consequence that the wave equation, Eq. (Q, becomes scale invariant which is very 
convenient from a numerical point of view since results for one characteristic length 
scale can easily be scaled to a different length scale. 
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For the dielectric material itself one likewise defines a material dispersion D m 
by solving the wave equation with e{f^_,uj) = s m {uS) — n 2 m (uj). From the simple 
homogeneous-space dispersion relation it readily follows that 

A d^n m 

D ™ = ~- C ^y- (6) 

Intuitively one might speculate that the two kinds of sources of dispersion simply add 
up and actually the approximation 

D « D w + D m (7) 

is used widely in the literature. While the approximation is useful in qualitatively 
understanding the zero-dispersion properties it is however also clear (see e.g. the work 
of Ferrando et al. |17|) that quantitative correct results requires either a self-consistent 
solution of the wave equation or some accurate perturbative method [181 119) . 

In this paper we use a fully self-consistent solution of the wave equation, Eq. Q . 
For the dielectric function we use the frequency-independent value e = 1 in the air- 
hole regions while we for silica employ the usual three-term Sellmeier polynomial 
description, 



(A) =n 2 m (A) = 1 + ^2^^ (8) 



where the absorption lines Xj and the corresponding strengths aj are given by 

Ai = 0.0684043 fjm, ai = 0.6961663, (9) 

A 2 = 0.1162414^m, a 2 = 0.4079426, (10) 

A 3 = 9.896161 fj,m, a 3 = 0.8974794. (11) 

Figure illustrates the typical dispersion properties of a photonic crystal fibre. 
The strongly negative material dispersion D m of silica below A ~ 1.3 tend to make 
the total dispersion D of standard fibres negative for A < 1.3 /im, simply because of the 
very weak waveguide contribution D w . However, photonic crystal fibres are contrary 
to this since the composite air-silica cladding is seen to provide the guided mode with 
a strongly positive waveguide dispersion D w which tends to shift the zero-dispersion 
wavelength Ao far below 1.3 /im towards the visible. While the material dispersion is 
fixed the waveguide dispersion varies strongly in the phase-space spanned by d and A 
and in this way the competition between waveguide and material dispersion becomes 
a powerful mechanism in engineering the zero-dispersion wavelength. 



3. A strand of silica in air — the ultimate limit? 



As mentioned in the introduction the zero-dispersion wavelength may be pushed to 
lower values by simply increasing the air hole diameter and decreasing the pitch. 
For d/A — > 1 this will to some extend effectively leave us with a single strand of silica 
surrounded by air. This limiting case has been emphasized previously in the literature 
and in Fig. we reproduce the zero-dispersion wavelength results reported by Knight 
et al. |15| . Generally, the strand of silica will have either two dispersion zeros or none, 
with the exception of the special case where it only supports a single dispersion zero. 
PCFs turn out to follow the same overall pattern and the existence of two dispersion 
zeros turns out to have interesting applications in super-continuum generation |20j . 
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The results in Fig. [3] leave promises for a zero-dispersion wavelength down to 
below 500 nm which will eventually also be the ultimate limit for silica based PCFs. 
However, as also indicated by the dashed line the large index contrast between air and 
silica in general prevents single-mode operation at the zero-dispersion wavelength. In 
the following we will study to which degree the photonic crystal cladding concept of 
PCFs can be used to circumvent this problem. 

4. Photonic crystal cladding as a modal sieve 

As demonstrated already by Birks et al. J3| the photonic crystal cladding of a PCF 
acts as modal sieve which may prevent localization of high-order modes to the core 
region. This so-called endlessly single-mode property has later been studied in great 
detail |141 1211 1221 and it was recently argued that the endlessly single-mode 
phenomena is a pure geometrical effect and that the PCF is endlessly single mode 
for d/A < 0.42 irrespectively of the fibre material refractive index The photonic 
crystal cladding thus serves to limit the number of guided modes and at the same 
time the guided modes will to some extend inherit the chromatic dispersion properties 
observed for the strand of silica in air ■ 

The problem of zero-dispersion wavelength versus pitch has previously been 
studied for PCFs with an infinite photonic crystal cladding PJ] demonstrating curves 
qualitatively resembling the curve in Fig. [3] Here, we extend that work to PCFs with 
a photonic crystal cladding of finite spatial extent. In particular, we study the effect 
of a varying number N of rings of air holes surrounding the core region. We also 
study the cut-off and leakage properties to explore the possibility for a single-mode 
PCF with a zero-dispersion wavelength in the visible. 

Our numerical solutions of the wave equation, Eq. QJ, are based on a finite- 
element approach which is described in detail in Ref. . For the calculation of the 
cut-off wavelength and the leakage loss we refer to Refs. [20 E) and references therein. 

The results of extensive numerical simulations are summarized in Fig. QJ First 
of all we notice that the zero-dispersion wavelength versus pitch has a curve-shape 
qualitatively resembling the result in Fig. [3] for a strand of silica in air. Furthermore, 
we see that the number N of rings of air holes has little influence on the zero-dispersion 
wavelength. In particular, the results for N = 10 are in full quantitative agreement 
with those reported in Ref. (23- On the other hand, the spatial extent N x A of the 
photonic crystal cladding is as expected seen to have a huge impact on the leakage 
loss |2E1 HZ] as seen from the red shading indicating the region with a leakage loss 
exceeding 0.1 dB/km. Furthermore, N has as expected little effect on the cut-off 
wavelength since the cut-off and the modal sieving is governed by the width (A — d) 
of the silica regions between the air holes [21] rather than the spatial extent N x A of 
the photonic crystal cladding. 

Finally, we note that by choosing d/A ~ 0.7 we may realize a PCF with a single 
zero-dispersion wavelength down to ~ 700 nm with the fibre being single-mode for 
longer wavelengths. We believe this to be the ultimate limit for silica-based PCFs 
having a photonic crystal cladding comprising a triangular arrangement of circular air 
holes. Such results have been demonstrated experimentally by e.g. Knight et al. |15j . 
In practice, the limit might be pushed slightly further toward the visible since real 
PCFs tend to have a slightly shorter cut-off wavelength compared to the expectations 
based on the ideal fibre structure [231 - Most likely, this tendency originates in the 
presence of scattering loss in real fibres which also acts in suppressing the high- 
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Figure 4. Zero-dispersion wavelength Ao versus pitch A for different values of the 
normalized air-hole diameter d/A. Panels a), b), and c) are for N = 6, 8, and 10 
rings of air holes, respectively. Regions with a leakage loss larger than 0.1 dB/km 
are indicated by red shading. Similarly, the multi-mode regime is indicated by 
blue shading. 
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order modes even though they are weakly guided by the photonic crystal cladding. 
In order to push the zero-dispersion wavelength further into the visible one would 
have to tolerate guidance of high-order modes or alternatively employ somewhat more 
complicated designs involving a varying air-hole diameter throughout the cladding . 

5. Conclusion 

In conclusion we have studied the zero-dispersion wavelength Ao in silica based 
photonic crystal fibres with special emphasis on the interplay with the cut-off 
wavelength and leakage loss. In the large parameter space spanned by the air-hole 
diameter d and the pitch A we have identified the values facilitating the shortest 
possible zero-dispersion wavelength (~ 700 nm) while the fibre is still single-mode for 
longer wavelengths. 

We believe that our Ao-maps are an important input for the efforts in designing 
nonlinear photonic crystal fibres with still shorter zero-dispersion wavelengths for 
super-continuum generation in the visible. 
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